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SUMMARY 

A i study of the effect of the frequency of the lowest 
winy structural mode on the airplane cent er-ef-y rarity 
dynamic-response factor was made by employing 
simplified transfer functions. It was found that the 
simplified transfer function adequately predicted the 
maximum value of the incremental normal-load- 
factor response at the airplane center of y rarity to 
isosceles triangle pulse elevator inputs. 

The results of the study are presented in the form of 
preliminary design charts which give a comparison 
between the dynamic-response factors of the semi- 
rigid case and the airplane longitudinal short-period 
case and between the dynamic-response factors of the 
semirigid case and the steady-state value of the air- 
plane longitudinal short-period response. These 
cha rts can be used to est imate the first-order effects of 
the addition of a wing-bending degree of freedom on 
the short-period dynamic-response factor and on the 
maximum dynamic-response factor when compared 
with the steady-state response of the system. The re- 
sults show that a structurally damped frequency 
greater than six times the short-period damped fre- 
quency will not affect the dynamic-response factor of 
the semirigid short-period response at the airplane 
center of gravity and that, when the frequencies are 
equal , the semirigid dynamic-response factor may be 
as much as 1.6 times that of the short period. The 
results also show that the maximum dynamic-response 
factor can be as much as 2./+ times the steady-state re- 
sponse of the system , depending upon the ratio of the 
natural frequencies of the structural and short- period 
modes and upon the damping of the two modes. 

INTRODUCTION 

As airplanes have increased in size, speed, and 


flexibility, analysis of the loads, stresses, and de- 
flections associated with the longitudinal shorl- 
period mode has become increasingly more com- 
plex. This complexity results from the need to in- 
elude not only the aeroelastic effects associated 
with structural deformation but also the dynamic 
effects of structural vibration. Considerable effort 
is currently being expended in the field of dynamic 
analysis and it lias become customary to express 
the dynamic effects of loot li aeroelasticily and 
st ructural vibration in terms of a dynamic-response 
factor which relates the dynamic response of the 
airplane to its steady-state response. The effects 
of flexibility are generally associated with a specific 
response at the center of gravity of the airplane, 
especially in the preliminary design stages; how- 
ever, these effects at other points on the airframe 
(such as a wing-tip deflect ion or a strain in a part ic- 
ular structural member) are often of interest. 

The present-day use of thin high-aspect -ratio 
wings on large high-speed airplanes has resulted in 
a lowering of the frequency of the wing structural 
vibratory modes. As a consequence of this re- 
duction in stiffness, the frequency of the lowest 
wing vibratory mode is approaching the frequency 
of the airplane short -period mode. The proximity 
of the frequencies of these two modes has a pro- 
nounced effect on the airplane dynamic-response 
factor. Although this effect lias been known quali- 
tatively for some time and studies of specific con- 
figurations have been made, there has been no 
simple' numerical guide' for estimating the effects 
of this design trend. Possibly, this lack is a natural 
consequence of the nature of the mathematical 
transfer functions which relate the airplane center- 
of-gravity response to an incremental change in 


1 Supersedes XACA Technical Note 4250 by Carl R. TTuss and Tames J. Doncpan, I95S. 
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elevator angle. These transfer functions are of a 
type which is usually regarded as being more 
adaptable to specific studies than to generalization. 

The purpose of this study was to determine 
whether the results obtained by using the complete 
transfer functions could also be obtained to a high 
degree of approximation with related but greatly 
simplified transfer functions and whether this sim- 
plification was of such a nature as to permit 
generalization of the results. The present paper 
illustrates the nature and validity of the simplifi- 
cation of the transfer function used and assesses as 
to botli magnitude and trends the effect of the 
proximity of the frequencies of the lowest wing 
structural mode and the airplane short-period 
mode on the airplane incremental normal load 
factor at the airplane center of gravity. The 
results are summarized in the form of design charts 
whiclq it is believed, will be of value in the pre- 
liminary design stages of an airplane. 

SYMBOLS 

A Z h generalized non dimensional mass- 

coupling term between Z and h 
degrees of freedom, a xh fpSc 
A hh generalized nondimensional mass 

term of flexible-wing mode be- 
tween elastic wing and h degree 
of freedom, a hh /pSc 

A 9Jt generalized nondimensional mass- 

coupling term between 0 and h 
degrees of freedom, a 9h jpSc 2 
a Zh generalized mass-coupling term be- 

tween Z and h degrees of freedom, 

2 J ^ ‘ m ' “U ( y } 1 ~ s'd* Cv) } 

sings 

a hh generalized mass term of flexible- 

wing mode between elastic* wing 
and h degree of freedom, 

2j o 1 Im'AWyW 

-2S’JMU(y) +rAMy)] 2 )<!yo, 

slugs 

a& h generalized mass-coupling term be- 

tween 6 and h degrees of freedom, 

2 j [— S'rMy)-\ r,Mv) 

-m'jj'W+S'JMyWv; 

slug-ft 


b wing span along elastic axis, ft 

C U| .. . C» dimensional transfer-function 
coefficients for semirigid 
case 


O' i f O' 2 , . . . C f 9 nondimensional transfer- 



function coefficients for 
semirigid rase 

Cf 

force coefficient due to elastic-wing 
d (‘flee tion, FJqS 

c m 

airplane pitching-moment, coeffi- 
cient about the center of gravity, 
M/qSc 

Cy 

airplane normal-force coefficient at 
the airplane center of gravitv, 
F„lqS 

c 

”~T 

d!~ 

c 

wing mean aerodynamic chord, ft 

e k 

kinetic energy, ft -lb 

E, 

potential energy, ft-lh 

F 

force, lb 

Fy 

airplane normal force at center of 
gravity, positive downward, lb 

My) 

spamvise bending-mode shape along 
wing elastic axis 

My) 

span wise twisting-mode shape about 
wing elastic axis per unit tip 
bending deflection, radians/ft 

u 

acceleration due to gravity, ft/sec 2 

ii 

wing-tip deflection, hfc, chords 

h 

wing-tip deflection of elastic axis 
due to bending, positive down- 
ward, ft 

r 

section moment of inertia, m'^x 2 , 

slug-ftVft 


K^K l} Ko } \ dimensional rigid transfer-function 

K 7 ,K^,Kq J coefficients 
K'iJy'iJy':, \ nondimensional rigid transfer-func- 
K' 8) K'<j j tion coefficients 
K y radius of gyration about T-axis, 

chords 

Ic r(‘duced angular frequency, oicjV 

l longitudinal distance from airplane 

center of gravity to wing elastic 
axis (function of span wise loca- 
tion), positive forward, ft 
M pitching moment about center of 

gravity, ft -lb 
m Ay m h m w mass, slugs 

m' w section mass, slugs/fl. 
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An 

incremental normal load factor at 
airplane center of gravity (posi- 
tive upward), g units 

<1 

dynamic pressure, lb/sq ft; also, 
pitching angular velocity, radi- 
ans/sec 

S 

wing pi an -form area, sq ft 

S' 

section mass moment about elastic 
axis, m'tP*, slug-ft/ft 

s 

Laplace transform variable 

T 

2 7 r 

period of oscillation, — , sec 

Ct) 


duration of triangular input, sec 

t 

time, sec 

V 

velocity, fps 

X 

longitudinal displacement, positive 
forward, ft 

V 

lateral or span wise displacement, ft 

/_ 

vertical displacement of airplane 
center of gravity, positive down- 
ward, ft 

z 

vertical wing deflection of elastic 
axis due to wing bending, posi- 
tive downward, ft 

<x 

angle of attack positive when wing 
leading edge is up, radians 

r«, Tn • ■ 

. r 5 dimensional quasi -steady trans- 
fer-function c oeffic i on t s 

r'i, r'», . 

. . r' 5 nondimen sional quasi -steady 
transfer-function coefficients 

T 

dynamic-response factor at air- 
plane center of gravity, 

(A^-dyn)mar/ static 


incremental elevator deflection, pos- 
itive when trailing edge is up, deg 

6 

angle of pitch about airplane center 
of gravity, positive nose up, 
radians 

A 

sweep angle of elastic wing, deg 

M 

nondimensiomd airplane mass, 
m A /pSc 

£ 

damping parameter, percent of criti- 
cal damping 

p 

mass density of air, slugs /cu ft 


angle of twist of airfoil in plane 
perpendicular to clastic axis, posi- 
tive when wing leading edge is 
up, radians 

w 

angular frequency, radians/sec 


Subscrip ts: 
.1 


d dampen! 

dyu dynamic 

f fuselage; structural 

h flexible-wing degree of freedom 

max maximum 

n natural 

o exposed wing 

sp short period 

sr semirigid 

w wing 

Z vertical degree of freedom 

6 pitching degree of freedom 

Dots are used to indicate differentiation with 
respect to time; for example, 0=ddjdt. The sub- 
scripts a, 6, A, h } <j, IT, and 8 indicate differentiation 
with respect to the subscripts; for example, 
C Na ^dC N /da. 

GENERAL CONSIDERATIONS 

In the preliminary design stage of an airplane, 
the designer can, with presently available meth- 
ods, estimate the longitudinal short-period dy- 
namic-response characteristics of the center of 
gravity of a given configuration for rigid and 
quasi-steady airframes. The rigid airframe is de- 
fined in this paper as a structure that does not 
deform or vibrate, the quasi-steady airframe ns 
one which can deform but not vibrate, and the 
semirigid airframe as one which can both deform 
and vibrate. The problem that the designer is 
faced with in this preliminary design stage is the 
effect of the airframe vibratory modes (particu- 
larly those of the wing since it is usually the most 
flexible) on the quasi-steady airframe longitudinal 
short-period dynamic response. The methods 
available for calculating these effects are usually 
rather complex or require information which would 
probably not be readily available at this stage of 
the design. The designer needs, therefore, some 
means of estimating these effects which are simple 
and are based on parameters which would be 
available. 

In this paper such means are presented in the 
form of preliminary design charts which can be 
used to estimate the effects of the proximity of 
the frequencies of the lowest wing structural mode 
and the airplane short -period mode (quasi-steady 
case) on the dynamic response at the center of 
gravity of the semirigid airplane. The design 
charts are based on the philosophy that in the 


airplane 
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preliminary design stage of a particular configura- 
tion the designer will be able to compute either 
the maximum incremental normal load factor for 
the quasi-steady case or the steady-state value of 
the incremental normal load factor of the quasi- 
steady case. The charts are restricted to estimat- 
ing the effects of only the lowest vibratory wing 
bending mode on the incremental normal-load- 
factor response at the center of gravity to elevator 
control inputs. All other structural parts are con- 
sidered rigid. The charts are further restricted to 
a comparison of dynamic-response factors which 
are defined as 

y~ (Aw ayn) static (1) 

whore (Ari dl , n ) max is the maximum amplitude of the 
first peak of the time history of the incremental 
normal load factor at the center of gravity and 
A n statu is the steady-state amplitude of the time 
history of incremental normal load factor. 

METHOD OF ANALYSIS 

The procedure followed in this paper for study- 
ing the effect of the proximity of the frequencies 
of the lowest structural wing mode and the short- 
period mode on the incremental normal -load 
dynamic-response factor at the airplane center of 
gravity was patterned after that of references 1 
and 2. Dynamic systems representing the incre- 
mental normal-load-factor response at the center 
of gravity to an elevator input and defined mathe- 
matically by transfer functions were excited by 
various isosceles triangular inputs and the maxi- 
mum values of the resulting time responses were 
expressed as ratios to the steady-state response 
factors. This procedure was followed for systems 
having the quasi-steady mode coupled with a 
structural mode (semirigid case) and for the quasi- 
steady mode alone for a wide range of configura- 
tions and frequencies and dampings of the two 
modes. The dynamic-response factors thus ob- 
tained for the system with two modes were then 
expressed as ratios to those obtained for the sys- 
tem with one mode to determine the effects in 
question. 

Although triangular inputs were used in this 
study, it is believed that comparable results would 
be obtained for other shapes of pulse-type inputs 
since the process of expressing t lie semirigid results 
as ratios to the short-period results tend to elimi- 


nate the effects of different -sli aped inputs. Isos- 
celes triangle inputs were used in this paper for 
the following reasons: they approximate in shape 
severe pilot-imposed inputs; their frequency con- 
tent could he easily varied by changing tlieir 
duration T t \ their frequency content could be 
made sufficient to excite the wing structural mode; 
and they could be easily handled mathematically 
both by manual calculation and by automatic elec- 
tronic calculation. 

For existing airplanes with lugh-aspect-ratio, 
thin, flexible wings, the lowest structural fre- 
quency is usually associated with wing bending 
and, therefore, wing bending was selected as the 
lowest structural mode for this investigation. The 
theoretical system chosen for this study consisted 
of three degrees of freedom: freedom in pitch, 
vertical translation, and wing bending. The equa- 
tions of motion developed by Lagrange’s method 
describing this system have been previously estab- 
lished and are presented in reference 3. For con- 
venience they are also restated in appendix A of 
this paper. 

The assumptions made in this study included 
the following: linearity, no change in airplane for- 
ward velocity, small perturbations, and rigidity of 
the fuselage and tail assemblies. Those assump- 
tions may be summarized by the assumption that 
the motions of an aircraft with flexible wings arc 
described by the equations given in appendix A. 
It was further assumed that the aircraft is stati- 
cally and dynamically stable longitudinally, that 
is, that the aircraft short-period mode and struc- 
tural mode are oscillatory and are damped. 

As mentioned earlier in this section, the dynamic 
systems used in this paper were defined mathe- 
matically by transfer functions relating the incre- 
mental normal load factor at the center of gravity 
to an incremental elevator input. Some of the 
terms in the transfer functions could be eliminated 
with small loss in accuracy and the analysis was 
made by using these simplified transfer functions. 
In order to show this relationship, it is first 
necessary to define the complete transfer functions 
and then demonstrate the simplifications that 
can be made to obtain the simplified but practical 
transfer functions. Hereafter in this paper the 
word "complete” will refer to transfer functions 
containing all the terms and the word “simplified” 
will refer to the transfer function with some of its 
terms omitted. 
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COMPLETE TRANSFER FUNCTIONS 

The complete transfer functions relating the 
incremental normal load factor at the center of 
gravity to the incremental elevator angle input 
for the semirigid case, the quasi-steadv case, and 
the rigid case were obtained from the equations 
of motion given in appendix A. 

Semirigid case. — The transfer function for the 
semirigid case which defines a system that has 
both wing quasi-steady deformation and wing 
vibration is, in nondimensional form, 

Ar* , , C'^+C'^+0' 7 s*+C'sS+C if » 

A5 e C ,s* 3 ~+ C V+ c\ 

where the definitions of the C f coefficients are 
given in appendix B. In dimensional form the 
transfer function may be written as 

An C^+C^+(W+WC, 

A5 e [S) - + C v s * 3 + CV+ C 3 *+ r 4 w 

where the conversion factors of C f to C arc given 
in appendix C. The static value of this transfer 
function is seen to be C ,JC A . The characteristic 
equation may be factored into two quadratic 
equations by GraefiVs method and written as 

An c 5 « 4 + r> 3 + cy + o+ r, 

aC S) [8*+2S 9p M'J 1+ M. P *)[#+2irMj*+ Mf 2 } 

( 4 ) 


Quasi-steady case. — As indicated in reference 4 
by letting rates of wing-tip deflection D 2 h=Dh ~ 0 
in the equations for the semirigid case, the transfer 
function for the quasi-steady case may be formed. 
In this case the wings can deform but do not 
vibrate. The transfer function for the quasi- 
steady case may be written as 


A n . . r 3# 2 h 4 s 5 

a^ w_ 5 2 +r 1 .s*+r 2 


( 5 ) 


where the F coefficients are defined in appendixes 
B and C. The static value of this function is 
seen to be r 5 /T 2 . It is interesting to note that 
the static value of the semirigid case is equal 
to the static value of the quasi-steady case 
CJC a = T 5 /T 2 , since r 0 r' 2 -C 0 C' 4 and r 0 r' 5 =f^V 
Rigid case. —By letting D 2 h=Dh—h — i) in the 
equation for the semirigid case, the transfer 
function for* the rigid case may be formed and 
written as 


5 

A // , , _ K 9 s 2 + A> -f K 7 ^ 

a£ w “ s'+Ks+ru w 

where the K coefficients are defined in appendixes 
B and O. Although the transfer functions of the 
rigid case and the quasi-steady case have the 
same form, the transfer coefficients of the two 
cases differ in that the transfer coefficients of the 
rigid case arc modified by the effects of w 7 ing 
deformation to give the quasi-steady-case transfer 
coefficients. The static value of this function 
is A 7 /A 2 . 


SIMPLIFIED TRANSFER FUNCTIONS 

Iii order to determine whether the number of 
terms in the complete transfer functions could be 
reduced, studies of 15 airplane configurations 
were made on an electronic analog computer by 
using the complete transfer functions for the 
semirigid case (eq. (4)) and the quasi-steady 
case (eq. (5)). The configurations used covered 
a range of wing sweep angles from 0° to 60°, of 
ratios of wing mass to airplane mass from 0.15 
to 0.50, of airplane center-of-gravity positions 
from 0.25c to 0.45c, and of dynamic pressures 
from 100 to 800 pounds per square foot. These 
studies indicated that some of the terms in the 
numerator of the transfer functions did not con- 
tribute appreciably to the maximum value of the 
time history of incremental normal load factor 
for triangular inputs but merely aided as phase 
shifters and thus were not required for the pur- 
poses of the present study. Typical results of 
these studies are shown for the semirigid case in 
figures 1 and 2 and for the quasi-st eady case in 
figures 3 and 4. 

Semirigid case.- — The contribution of the nu- 
merator terms of the semirigid transfer function 
with each numerator coefficient equal to unity is 
shown in figure 1(a). In figure 1(b) the con- 
tribution of each of these same numerator terms 
is shown for typical values of the coefficients. 
From plots such as these it is seen that the C 9 term 
makes the most important contribution to the 
maximum value of the incremental normal load 
factor. 

Calculations of incremental normal-load-factor 
time response to isosceles triangle inputs were 
then made by using only the C 9 term in the 
semirigid transfer function. These time histories 
were compared with time histories’ obtained from 
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Time, t, sec Time . '• sec 


(a) Contribution of numerator terms with coefficients 

equal to unity, 

(b) Contribution of numerator terms for typical values of 

the coefficients: (79=563.25; (7s — 50.376; (77=0,604; 
C 6 = 0.504; and C b =- 0.0394. 

Ficukk 1. Contribution of individual numerator terms 
uf the semirigid transfer function 

Aw . v 4 CqS* 4 C7$ 8 4 C&s -f C o 

A5~ [S) ~ « 4 + C^+C 2 s 2 +C 3 «+C 4 “ 

to the complete response to a unit-amplitude isosceles 
triangle input with 77=1.0, C|=15.4414, (72=116.8380, 
= 350.6639, and ^=554.5209. 


(a) Mo 'Vs apart; (w d ) / /(aij), p =2.44. 

(b) Modes together; («*),/(«*),„- 1.0. 

Ficttre 2. Comparison of the response to a unit-ampli- 
tude isosceles triangle input with 77 1.0 obtained from 

using the simplified semirigid transfer function 



C» 

* 4 - fC,^ + C^+C 3 s+C 4 


and the* complete semirigid transfer function 


Aw (V 4 C 6 .s 3 4 CV* 2 + C*s f C v 
AS e W a* f CiS* 4 C 2 S 2 + C,s 4 C 4 ‘ ’ 


tlie complete transfer function (eq. (4)), to de- 
termine how well the simplified transfer function 
((7 9 term only in the numerator) described the 
maximum value of the time response of incre- 
mental normal load factor for triangular inputs. 
Typical comparisons are shown in figure 2(a) for 
the case when the frequencies of the modes are 
different and in figure 2(b) for the case when 
the frequencies of the modes are equal. 

On the basis of such computations it was 
determined that the complete semirigid transfer 
function (eq. (4)) could be reduced to 


( w Jsp 2 ] [* 2 + (w rt )/S -f (w«)/ 2 ] 

( 7 ) 


and still adequately describe the maximum value 
of the time history of incremental normal load 
factor for triangular inputs. The use of the 
word “adequately” in this paper means generally 
to within about 3 percent and rarely more than 
about 10 to 15 percent. 

Quasi-Steady case. — A similar procedure was 
used to determine the contribution of the terms in 
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the numerator of the quasi-steady transfer func- 
tion (cq. (5)) to the maximum value of the time 
response of the incremental normal load factor. 
In figure 3(a) the contribution of the numerator 
terms of the quasi-steady transfer function with 
each numerator coefficient- equal to unity is 
shown. In figure 3(1)) the contribution of each 
of these same numerator terms is shown for 
typical values of the coefficients. In this case it 
is seen that r 5 is the important term. Typical 
comparison of a time history obtained from (lie 
reduced transfer function (r 5 term only in the 
numerator) with that obtained from the complete 
transfer function (eq. (5)) is shown in figure 4. 



(a) Contribution of numerator terms with coefficients 

equal to unity. 

(b) Contribution of numerator terms for typical values 

of the coefficients: r,- 7.093; Pp— 0.054; and 

r 3 = -0.0355. 

Ficttre 3.— Contribution of individual numerator terms 
of the quasi-steady transfer function 

Sn , . r :i s 2 + r 4 .s-r Ps 

to the complete response to a unit-amplitude isosceles 
triangle input with T,— 1.0, 3.4328, and P 2 =G.993f. 

49S505- 50 - -2 


From comparisons such as that shown in figure 
4, it was determined that the complete quasi- 
steady transfer function (oq. (5)) could bo re- 
duced to 

and still adequately describe the maximum value 
of the time history of incremental normal load 
factor for triangular inputs. 

Semirigid short-period case. -Since the denomi- 
nator of equat ion (8) does not equal the short - 
period part of the denominator of equation (7), 
it was found convenient to define another transfer 
function. This transfer function will be called 
the semirigid short-period case and is defined as 

A n ( s ) = — (q\ 


where the denominator of equation (9) is identical 
to the short-period part of the denominator of 

C 

equation (7) and the constant A is equal to . 2 » 

The use of this semirigid short-period transfer 
function as a basis of comparison rather than the 
quasi-steady case reduced the computations to 
practical proportions. If the quasi-steady ease 
had been used, it would have been necessary to 
estimate a new set of derivatives which make up 



Fioure 4. — Comparison of the response to a unit-ampli- 
tude isosceles triangle input with 7*= 1.0 obtained from 
using the simplified quasi-steady transfer function 


-IL = £» 

and the complete quasi-steady transfer function 

An = ] r 8 g» + r 4 «+ r a 

M* s 2 -I-riS+fa ' 
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the transfer-function coefficients for each new 
configuration and flight condition. However, by 
defining the semirigid short period, it was neces- 
sary to choose only the damping and frequency 
of the t wo modes without regard to the derivatives 
which determine these parameters. 

Actually, the semirigid short-period case is 
practically equal to the quasi-steady case since t lie 
damping and frequency of the two cases are almost 
the same for a wide range of configurations and 
q values. (See figs. 5 and 6.) A comparison 
of the natural frequencies of the semirigid short- 
period case with those of the quasi-steady case for 
a wide range of configurations and q values is 
shown in figure 5. A similar comparison of the 
damping of the semirigid short-period case with 



Figure 5.- Comparison of the* natural frequency of the 
semirigid short-period mode with that of the quasi- 
steady mode at a conler-of-gravity location of 0.25c. 
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Figure 6. — Comparison of critical damping of the semi- 
rigid short-period mode with that of the quasi-steady 
mode at a centcr-of-gravity location of 0.25c. 

that of the quasi-steady case is shown in figure G. 
The points shown in figures 5 and 6 were computed 
from the data of reference 3. The data of figure 6 
indicate that at the higher values of dynamic 
pressure the damping of the short-period case is 
greater than that of the quasi-steady case. There- 
fore the maximum value of the response as well 
as t lie maximum dynamic-response factor obtained 
from the short-period ease would always he less 
than that of the quasi-steady case. Tlius the 
ratios of maximum dynamic-response factors ob- 
tained by comparing the semirigid case with the 
short-period case would always he greater than 
(on the conservative side) or the same as those 
obtained by comparing the semirigid case with 
the quasi -steady case. 

Rigid case. -The complete rigid transfer func- 
tion (eq. (6)) could also be reduced to 



200 400 600 800 1,000 

Dynamic pressure, q y Ib/sq ft 
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(A 

A 8, w s*+K v s+K 2 


( 10 ) 


in order to define the maximum value of the time 
history of incremental normal load factor for tri- 
angular inputs. 


TYPICAL CALCULATIONS ILLUSTRATING METHOD 

The method used in this study and the prepara- 
tion of the desired preliminary design charts in- 
volved a large number of calculations and plots, 
typical samples of winch are shown in figures 7 
and 8. 

In order to obtain the maximum possible dy- 
namic-response factor for the range of the vari- 
ables, it was first necessary to calculate the time 
response to triangular inputs of varying duration 
(different frequency content) for each system (a 
particular combination of the variables). A 


sample of these calculations is shown in figure 7 
for both the semirigid short-period and the semi- 
rigid cases. Some of these computations wore 
carried out on automatic electronic computing 
equipment, some on desk-type computers, and 
some were carried out by using the tables of refer- 
ences 5 and G in conjunction with automatic 
electronic computing equipment. The dynamic- 
response faetoi- defined previously as equation (1) 
w T as determined for each ease by picking the value 
of the first peak of the time histories (see, for 
example, fig. 7 for ^=0.4) and dividing it by the 
static value for the particular system being con- 
sidered. These results were then plotted against 
the period ratio T f /T sp (ratio of the time base of 
the input to the natural period of the short-period 
mode) in order to determine the maximum dy- 
namic-response factor for each case. A typical 
plot of this procedure is shown in figure 8. The 



Figure 7. — Typical center-of-gravity incremental normal-load-factor time responses of the semirigid short-period and 
semirigid cases with a frequency ratio (to »),/ («„) of 1.9 to isosceles triangle inputs. («,) =3.162: £ =0.38: 

W/-0.0; ^ -0.02. 
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Period ratio, T t /T sp 

Fiourb 8. — Typical plot of dynamic- response factor y 
against period ratio TJT, P for the semirigid short- 
period case and semirigid case with a frequency ratio 
W;/W iP of 1.0. W <p = 3.162; 0.38; («*),= 6.0; 

€/-■ 0 . 02 . 

data of figure 8 are for the same cases as those of 
figure 7, the four points shown in figure 8 having 
been computed from the results shown in figure 7. 

From plots such as that shown in figure 8, if 
was possible to ascertain the magnitude and trends 
of the effects of the proximity of the frequencies 
of the lowest structural mode and the airplane 
longitudinal short-period mode on the incremental 
normal-load dynamic-response factor at the air- 
plane center of gravity. The plotting of these 
calculations resulted in the desired preliminary 
design charts. 

RANGE OF VARIABLES 

The results of this study are believed to be valid 
over a range of variables as follows: dynamic 
pressure from 100 to 800 pounds per square foot, 
wing sweep angles from 0° to 60°, ratios of wing 
mass to airplane mass of 0.15 to 0.50; centor-of- 
gravity location from 0,25 to 0.45 moan aoio- 
dynamic chord, ratios of damped wing lowest 
structural frequency to damped airplane' longi- 
tudinal short-period frequency from 1 to 15, and 
damping of the wing lowest structural mode and 
airplane longitudinal short-period mode from 0 to 
95 percent of critical damping. 

RESULTS AND DISCUSSION 

The results of this paper are summarized in the 
form of preliminary design charts. As mentioned 
previously, those design charts wore obtained fiom 
plots such as those of figure 8 covering a complete 
range of combinations of the variables. 


The chart given as figure 9 is a plot of the 
ratios of maximum dynamic-response factors 
(y,f)maxKy,p)max against tho ratio of the slruoUiral- 
modc' natural frequency divided by the semirigid 
short-period natural frequency, the damping of the 
structural mode being held constant at 2 percent 
of critical damping. The data of figure 9 were 
reduced to a more compact form by converting 
the abscissa to a ratio of the damped frequency of 
the structural mode and the damped semirigid 
short-period frequency. This simplification is 
given as the design chart shown in figure 10. 

The design chart shown in figure 10 indicates 
that, if the damped structural frequency and the 
damped semirigid short-period frequency are 
equal, the maximum dynamic-response factor of 
the semirigid case will he about 1 .6 times the value 
of the maximum dynamic-response factor for the 
semirigid short-period case. It can also be seen 
that, when the ratio of the damped structural fre- 
quency and the damped semirigid short-period 
frequency is greater than about 0, there is no in- 
crease in the maximum dynamic-response factor 
of the semirigid case over the value for the semi- 
rigid short -period ease. 

Increasing the damping of the short -period mode 
of the semirigid case while holding the damping of 
t J io structural mode constant also results in a de- 
crease in the maximum dynamic-response factor 
of the semirigid case. This effect can be seen from 
the results shown in figure 9. 

In order to investigate the effect of structural 
damping on the airplane dynamic-response factor, 
calculations were made for a semirigid short- 
period mode with a natural frequency of 3.1G2 
radians per second and a damping of 38 percent 
of critical damping coupled to a structural mode 
having variable damping of 0 to 95 percent of 
critical damping and a damped frequency equal 
to the semirigid short -period damped frequency 
(2.920 radians per second) and equal to 5 times the 
semirigid short -period damped frequency. Py- 
namic-response-f actor ratios for these eases are 
plotted against critical damping of the structural 
mode in figure 11, The result shown in figure 11 
indicates that, for a given value of damping of the 
short -period part of the semirigid case, an increase 
in the damping of the structural mode results in a 
decrease- in the maximum dynamic-response factor 
of the semirigid case. 
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FintTHK 9.- —Effect- on the ratio of maximum dynamic-response factors of the proximity of the structural natural fre- 
quency to the semirigid short-period natural frequency. ^—0.02. 


Titus, the data of figures 9, 10, and 11 indicate 
that, for a given frequency of the short-period 
mode of the semirigid case, an increase in the fre- 
quency and/or an increase in damping of the struc- 
tural mode and/or flu increase in damping of the 
semirigid short-period mode result in a decrease of 
the maximum dynamic-response factor of the 
semirigid case. Figures 9, 10, and 11 would be 
the ones used by a , designer in order to obtain an 
estimate of the increase in the incremental nor- 
mal-load short-period dynamic-response factor at 
the airplane center of gravity due to the proximity 
of the frequency of tine lowest wing structural 
mode to that of the airplane longitudinal short- 
period mode. Use of these figures presumes, as 
mentioned earlier, that the designer would be able 
to estimate the maximum longitudinal short- 
period response and would have an estimate of the 


lowest wing structural frequency and damping of 
the lowest wing struct ural mode. 

Another design figure which may be useful is 
one which gives the effect of the proximity of the 
structural natural frequency to the short-period 
natural frequency on the maximum semirigid 
dynamic response when compared with the semi- 
rigid short-period static value. This result was 
easily obtained by plotting the semirigid maximum 
dynamic-response factor for each case (obtained 
from plots such as fig. 8 and noting that, as pointed 
out earlier, the static value of the semirigid and 
semirigid short-period cases are equal) against the 
ratio of the structural natural frequency to the 
semirigid short-period natural frequency. Such 
a plot is presented as figure 12. 

The designer could use the chart given in figure 
12 under the same restrictions as wore mentioned 
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Figure 11, — Effect of damping of the structural mode on 
the ratio of maximum dynamic-response factors for two 
ratios of the damped structural and semirigid short- 
period frequencies for semirigid short-period damping 
of 0.38. 

for the previous charts with one exception. This 
exception is that in using this chart the designer 
would need to know only the airplane longitudinal 


short-period steady-stale response rather than the 
maximum short-period response. 

The effect of the input-time base on the dy- 
namic-rcsponsc-factor ratio may also be of interest 
and can be determined from plots such as that of 
figure 8. In this case, rather than express the, 
maximum values of the dynamic-response factor 
as ratios, the values of the dynamic -response fac- 
tor of the semirigid and short-period frequencies 
are expressed as ratios at specific values of the 
period ratio TijT sp and are plotted against the 
period ratio. Typical plots of this dynamic- 
response-factor ratio are shown in figure 13 for 
three values of short-period damping. 

Tlie base of the input that gives the maximum 
dynamic-response factor is, of course, different for 
each case, depending on the damping of the two 
modes. It was usually greater than about 0.7 
of the natural period of the short period for all the 
cases studied in this paper. Examination of plots, 
such as those shown in figure 13, indicate that, 
when compared for the same triangle base, the 


EFFECT OF WING BENDING ON AIRPLANE DYNAMIC-RESPONSE FACTOR 


13 



Figure 12. Effect of the proximity of the structural natural frequency to the semirigid short-period natural frequency 

on the semirigid maximum dynamic-response factor. £ f =0.02. 


highest ratio of dynamic-response factor for a fre- 
quency ratio of 1.0 will he obtained from triangu- 
lar inputs with a base equal to 0.G to 0.8 of the 
natural period of the short period. For frequency 
ratios greater than 1.0, the ratio of dynamic- 
response factors is greatest for triangles with a 
base equal to less than 0.1 of the natural period of 
the short period. Thus, it is difficult to pinpoint 
a specific triangle base as being the one giving the 
most severe results. 

Since airplanes operate at flight conditions (alti- 
tude, airspeed, center-of-gravity location) which 
are constantly changing, the frequency ratio for a 


particular configuration will not be constant. 
Present-day large high-speed airplanes with thin, 
high -aspect-ratio, flexible wings are operating in 
the frequency-ratio range of roughly 4 to 10. 
The conditions for winch the frequency ratio will 
he a minimum depends somewhat on the configu- 
ration hut, in general, operations at low altitude, 
high airspeed, and forward center-of-gravity posi- 
tion should result in the lowest frequency ratio. 
This effect can he seen in figure 14 where the effect 
of dynamic pressure and airplane configuration 
on the proximity of the damped frequency of the 
structural mode to that of the short-period mode 
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(c) *.,-0.05; V=0.02. 

FiorRK 13. — Typical plots of the effect of tile period ratio 
oh tin- ratio of dynamic-response factors for various 
frequency ratios, 

is given. The data of this figure wore converted 
from the data of reference 3. 

The data of figure 14 show that, for airplanes 
with unswept wings or wings with very little sweep, 
the frequencies of the modes are brought into 
closer proximity by an increase in the dynamic 
pressure or, for a given dynamic pressure, by 
moving the center of gravity forward. For 
wings with sweep angles greater than about 30°, 
these trends of the effects of dynamic pressure 
and centcr-of-gravity position on the proximity 
of the frequencies of the modes are the same. 
For these eases, however, the aeroelaslic effects 
caused by increasing the dynamic pressure usually 
cause the short-period mode to become statically 



Fjocrk 14 . — Effect of dynamic pressure q on the frequency 
ratio (u d ) ,/(«<*) iP tor various wing swoop angles, wing- 
mass-to-airplane-mass ratios, and center-of-gravity 
locations. 

unstable (indicated in fig. 14 by the frequency 
ratio going to infinity) betore the frequencies 
of the modes can be brought together. Thus, 
for a given configuration the operating conditions 
will determine the relative proximity of the 
frequencies of the two modes and at what point 
on the abscissas of the design charts the airplane 
is operating. 

It is well to emphasize that the preliminary 
design charts given arc only meant to give first- 
order effects and to apply only to systems which 
are statically and dynamically stable. Further- 
more, since the curve given in figure 10 is an 
envelope of the maximum values of the converted 
data of figure 9, it will normally give conservative 
values of the ratio of maximum dynamic-response 
factors. Finally, for a particular design problem 
a detailed analysis including all the variables 
should he made if the “rule-of-thumb” value for 
the ratio of maximum dynamic-response factors 
given by the chart indicates the possibility of a 
dangerous situation. 
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Figure 14. — Concluded. 

CONCLUDING REMARKS 

The results of this study of the effect of the 
frequency of the first wing bending mode on the 
airplane dynamic-response factor indicated that 
the maximum conter-of-gravity load-factor re- 
sponse to a triangular-shaped pulse elevator input 
could be adequately determined by using a 
simplified transfer function for the semirigid and 
quasi-static cases. The use of the short-period 
part of the semirigid transfer function as a basis 
of comparison gave results which were either equal 
to or on the conservative side of those that would 


have been obtained from the quasi-steady transfer 
function. 

As a result of the reduction in t lie number of 
terms obtained by using the simplified transfer 
functions, it was possible to construct design 
charts which provide trends and rule-of-thumb 
estimates of the effect of the frequency of the 
first wing-bending mode on the airplane dynamic- 
response factor. The charts show that the maxi- 
mum dynamic-response factor for the semirigid 
case will be l.fi times that of the short-period case 
when the damped frequencies of the structural 
mode and short-period mode are equal. Further- 
more, when the frequency ratio is greater than 
about C>, a lightly damped structural mode has 
little or no effect on the dynamic-response- 
factor ratio at the airplane center of gravity, and, 
as the damping of the structural mode increases, 
the frequency ratio at which the structural mode 
has negligible effect also decreases. Finally, the 
charts indicate that the semirigid maximum 
dynamic-response factor can be as much as 2.4 
times the steady-state value of the system, 
depending on the damping of the structural and 
short-period modes and on the ratio of the natural 
frequencies of the two modes. 

The dynamic-response factor for a particular 
configuration will vary with the operating condi- 
tions (principally with dynamic pressure) but 
should be a maximum at flight conditions of low 
altitude, high airspeed, and forward eentcr-of- 
gravily position. 

It should be repeated that for a particular 
design problem a detailed analysis should be 
made if the rule-of-thumb value given by (lie 
design charts indicates the possibility of a dan- 
gerous situation. 

Langt-et Research Center, 

National Aeronautics and Space Administration, 

Langley Field. Va.. February 21 , 1058. 



APPENDIX A 

EQUATIONS OF MOTION 


The equations of motion used in this paper will 
be given here briefly for the convenience of the 
reader. A complete development of these equa- 
tions is shown in references 3 and 4. The equat ions 
are derived on the basis of the Ijagrangian equa- 
tion: 


(I /dEA dE k dE p _ F 
dt\bQJ 


(Al) 


where 


E k kinetic energy 
E p potential energy 
Q generalized coordinate 


F q generalized force 

The three generalized coordinates used are: 

Z vertical translation 
6 pitching velocity 

h displacement of wing t ip due to bending of 
elastic wing 

For an unswept wing the flexible-wing mode 
shape consists of bending /,(}/); and for a swept 
wing the flexible- wing mode shape consists of 
bending /,(?/) combined with twisting per unit 
bending deflection at the wing tip /*(?/). The 
span wise bending is usually assumed to be para- 
bolic and the spanwisc twist, linear. 


The non dimensional equations of motion thus derived are. 


2 M 7?(«-0) +2 A zh mT-C s a-T)a 

2vK r t m+2A lh mi-r m a-naQn mi> ^-r)e(^ 2 (\ n ^-H(C mH )-nii(C ni>H )=c mM 'A& e 

2 A M ,mT\ 2 A Xh D(a-e) + 2 A H m+ 2 Ajm- C r a-D9 (’ cJ-H(C, h ) -DTI (CV M ) - C r „M t 

and by definition 

\n V /DO Da\ 

A5, fl \A6, A5 C / 


(A2) 


(A3) 


(A4) 


(A5) 


Simultaneous solution of equations (A2), (A3), (A4), and (Ao) results in the semirigid transfer function 
given as equation (3). 

1G 


APPENDIX B 

DEFINITION OF TRANSFER-FUNCTION COEFFICIENTS 


Tlu* transfer coefficients used for the analysis are defined in this appendix. 


SEMIRIGID CASE 


The coefficients for the semirigid ease arc 


r^A(„Kyh\ m -A^) (2^ ev, h .)-2(,l M .l z/i C m;ja + 

r " ,= A{ 4/jri '*^"( rv - + 2 . l M C „, j (2 M -Ic!v- ,„) 

-\-(4At h A zh -\-A hh O mi) ^ 2/i— g f. 'v, J ^+QuA^m /)a +2/l« A ^2^— - 

+ \.\ e/ ,(A ei /\- A xtl r m j +A ih c mD c IVl>H +2A^r mt y 

{wcvi.. (2,-1 r. m ) w»jtw„-r,.r % 4 2.i„r,„ )+ Q r. c,„„ 

— 2nK y C Fh -\- 2A )h C' m ^j (~f x ~2 1A th C ff ' — - 0.\- I)H C mu ^ ^2A Z „ fV a ^ 

+(2.i M r,„ o -2A w r',. i -2ri z/1 a„ n „-2 r’ Vf )+e.jA„e^-A x ,(^ 

^ Cv^j^J+s.Ija (e % r. Vjw - ev a C „ UJf +2 f 2 ^ 2 . 1 ^+^ cv,)} 

^ 3= r;^ — ' 4 ^Ai- 2 ev„— ^ (-V- 2 '^)+r m/to (-2 M -- r jV> ^J .i m p i- 2 m av( e V( /V„ 

~~ c + 2 ^0 /i(O fla ev„— + 1^2 (, '«„(^a x/i -\-^ < v^ 2 M — ^ e Vi) „^ 

— (y "'/-Vn/f + 2 flm »f S H~ C*<P’*i>ii) ( 2/lz/i +2 ^ 2 e*„A ZA — Cmf'v u „ 

-2 'VrAi) (-2^-2 ^0+(° A„ .CV^-24*^) i 

^^=^{[2 ^ a ^-i C„„(- 2 m -2 rv.)] ( 2 /UP-r, // )+(^ a ^„ 

-<Vk») ( 2 - 1^+2 rv„)+[r,/'„,„ (-2 m-J ev,)]+2 r > r -,,ev„} 
e "«=^ (4A M *C #l -4A«^ M r., -4 M AV^ev s ) 

C ' t= gf’„ ^** f ‘ v * ( r "» B +C I , f )-A w r mj (Cv /to +CvJ+2/l M r Vj (2A Zh —C mDH )+2 (A z/I r, „ t} 

—■‘'U»Cv,)(2A zt +2 r^+ 2f„, s (A^Cy nH —2A z A) +2nK i }C Fl)H C K ^ 
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T— 4,/Vr*Cvj4j W ^+24 z *C MJ (Cv /W l £>■„) l~2vl /,/,(C,„ a Cv 4 CxJ m ^-\-2A dn (Cs H ^m i (>J Vj) 

<J ( o L 

4 ™ J )(^24 //l +2 ^'(,^ + 2 < - F nn^'"i (Cv n< 4- CvJ 

— rj Dili’s (^rn/ >a 4 ^J4-2< r \’j(M^iVV // — 4 i)/ ,r , ,„ w )J 

C’'s=-7r -f - 1/,/4 2 Cv 4 (^ m „„+ ^ r« v ) — -' l*4.-'-V',„ s (r v „ <( + r,,) + f mj p N<Pri>lI ~trS'xi>H 

<J ( U t. L - 

+ 2 ^(^4 ( '’v,)j4 Cv 4 (hV ( /4, w ( F “^ n 'n<i\ ^ ra <r)J 

4"2-lz* (Cx H C mi - C mii C N j) + (f4,,„Cv 4 — Cv„ 4 ) ( 2 A Zh -\- ^C F<i ) J- 

r" 9 =-p [ {^m a ^ K&— Cv„C„ 4 ) (-- *"' / '4 '" V,,) d ^/» 4 Cv„)] 

in o 

QUA ST-STEADY CASE 

Tlu> coefficients used in the analysis for the quasi-steady case are as follows: 

T 0 = 2nKy- (2^ — i r v „P (£A llh k--C r ^+A en C mn J^N H + *vKY 2 ^?M r N H + ‘ZA ei r mH {^IX — g r * D a) 

r'i=-p-(p- 2juAVCv a — 2 f’„ 4 ( 2^-2 Cv ,,„)+2 t7 *ua( — 2 ^~ 2 <0 ' v *)] C 2 ' 1 ""^ - ^ A ’") 

+24 6A (C» a Ov w C Na ( 2 A zh C mH ( 2 ,— 2 Cv^ — ^ ( 2 ^ 2 /J< \) P' ^ 5 ) 

-Cv„(2M^,- 2 r Fa 4 ^0„,)-2 (-4^+2 G 0} 

r% =I \ u {[> rv a r m ,+r„, a (- 2,-1 r*,)] <aAJfc*-c*)+ (r* r^-a/vj ( 2 A„+± r , t ) 

KV/™,, (-2,-2 Cv t )+2 'V '/*//} 

r r 3 =“, [ — -nfcy 1 c s 6 (2 ^ U/,fc 2 fV 7/ ) ~h 2^U/,( h ^ n h Cn/ »»//)] 

r ^^« + ^“2 r W(^ /te 4f >,)] (2,U--r F „) 

+24z fl (a j r V/ - Cv j r',„ ;/ )+(2.u A +2 '•'.) (^v/W-^/’vJ } 

r' 5 =— r [(Cv/V C m{ Cv a ) (24 mP— C'r H ) + C Fo (Cv^ mtf — C,„ s Cv„)J 
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RIGID CASE 

The transfer functions used in the analysis for the rigid case are as follows: 


1 




) 


K« — 2/i7 ?Y 2 ^2 

B r «A-+ r - (-‘'-I '*.)] 

K '^W, (-f-.C-r KV’„.) 

A "'=^ [ C ’"‘ ( 2 "-2 C *-)+ C - (- 2 --2 '“.)+2 «+'’-,)] 


K '*-yiZo (~ 2 ^' 2< Ai) 


APPENDIX C 

CONVERSION FACTORS FOR DIMENSIONALIZING THE NONDTMENSIONAL TRANSFER-FUNCTION 

COEFFICIENTS 


The conversion factors for dimensionalizing the 
nondimensional transfer functions for the semi- 
rigid case are as follows: 



Tlie conversion factors for dimensionalizing the 
nondimensional transfer functions for the quasi- 
steady case are as follows: 



The conversion factors for dimensionalizing the 
nondimensional transfer functions for the rigid 


case arc as follows: 


II 

* GO- 

K * (!•>■> 

*• G> 


K„=K ' 9 
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